Superfocusing by Nano-Shells 
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Recently Merlin and co-workers demonstrated, both theoretically and in the 
microwave range experimentally, subwavelength focusing of evanescent waves 
by patterned plates. The present paper extends these ideas and the design 
procedure to scatterers of arbitrary shapes and to the optical range of wave- 
lengths. The analytical study is supported by numerical results. The most 
intriguing feature of the proposed design is that, in the framework of classical 
electrodynamics of continuous media, focusing can in principle be arbitrarily 
sharp, subject to the constraints of fabrication. 

OCIS numbers: 310.6628, 160.4236, 350.4238, 240.6680, 250.5403, 050.1970. 

© 2009 Optical Society of America 
OCIS codes: 



5-H 



c3 



i 

-a 
a 
o 

o 



> 

(N 
(N 



oo 
o 



X 



Subwavelength focusing that circumvents the usual 
diffraction limit in optics has been a very active area of re- 
search, with a multitude of approaches explored in the lit- 
erature: negative-index lenses and guides, plasmonic par- 
ticles and cascades, superoscillations, time-reversal tech- 
niques and others ([l]-[6] and references therein). Re- 
cently the Merlin and Grbic groups [7, 8] (see also [9]) 
showed, both theoretically and in the microwave range 
experimentally, that patterned (grating-like) plates pro- 
duce subwavelength focusing of evanescent waves if the 
pattern contains significantly different spatial scales of 
variation. Conceptually, these patterns are related to 
surface profiles of near-field optical holography [10] . 

This Letter shows that the ideas of near-field focusing 
can be extended to patterned nano-shells and, further, to 
scatterers of arbitrary shape. The most intriguing feature 
is "superfocusing" - i.e. focusing that can in principle be 
arbitrarily sharp and strong, subject to the constraints 
of fabrication and availability of materials with desired 
values of the dielectric permittivity e. (It is also tacitly 
assumed that the size of the system is sufficiently large 
for electrodynamics of continuous media to be applicable; 
see e.g. [6] and references therein.) 

Let us suppose that an incident plane wave with a fre- 
quency to impinges on a scatterer coated with a plas- 
monic and/or dielectric layer. In general, the shape of 
the scatterer may vary; more importantly, the thickness 
and dielectric permittivity of the coating may be chosen 
judiciously, with the ultimate goal of nano-focusing the 
wave at a given spot. In practical applications, the role 
of the scatterer can be played e.g. by the apex of an 
optical tip, the active part of an optical sensor or by a 
nano-antenna. 

Even though the ideas are general and could be applied 
in 3D electro dynamic analysis and design, let us consider 
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Fig. 1. A thin shell (coating) can be capable of focusing 
light to an arbitrarily narrow spot if the angular variation 
of its dielectric function and/or thickness are judiciously 
chosen. See text for further details. 



for maximum simplicity the 2D case of a cylindrical scat- 
terer (particle) of radius r cy \ (Fig. 1). The usual complex 
phasor convention with the exp(— hot) factor is adopted. 

Several ways of formulating this electromagnetic field 
problem are available. First, one needs to distinguish full 
electrodynamic analysis vs. quasi-static (QS) approxima- 
tions valid for dimensions much smaller than the wave- 
length. The QS treatment has limitations, as it does not 
account for wave effects and short-range surface plasmon 
modes [11]. However, to fix ideas and for the sake of 
algebraic simplicity this Letter adopts the QS approx- 
imation. Full wave analysis is analogous but involves 
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Bcsscl/Hankel functions and their derivatives instead of 
polynomials in r (see expressions below). 

Second, the electrodynamic problem may be formu- 
lated in terms of the electric or magnetic field (conve- 
nient for s- and p-polarizations, respectively). Alterna- 
tively, two dual formulations are also available under the 
QS approximation: via the electric scalar potential and 
via the stream function \&, V x $ = eE see [12, 13] 
for some interesting implications of this duality. 

Let us endeavor to achieve sharp focusing with respect 
to the polar angle ^ at a certain radius rf > r cy \. For 
p-polarization (the two-component electric field in the 
cross-sectional plane and the one-component magnetic 
field H directed along the axis of the particle), let 



H out (r f ,<f>) = H f g((p-<p f ) 



(1) 



where H out is the magnetic field outside the shell, Hf is 
some amplitude at the focus, <f>f is a given angle and g(<f>) 
is a desired distribution of the field. For example, choos- 
ing g{(f>) or, alternatively, g'{<f>) as a sharp Gaussian peak 
will result in the corresponding peak in the magnetic field 
or the radial component of the electric field, respectively. 

This local behavior of the field can be extended, by 
analytical continuation, to the whole region outside the 
scattering shell. One way of doing so is by expanding 
the field into cylindrical harmonics whose coefficients are 
found from the Fourier transform of g. This analytical 
continuation will then determine the field distribution on 
the surface of the particle, and one needs to find the pa- 
rameters of the coating that would produce such a surface 
distribution. 

More specifically, in the QS limit the field inside the 
shell can be expanded into cylindrical harmonics: 

oo 

H in (r,<f)) = a " rH exp(in^), r < r cyl (2) 

n— — oo 

where a n are some coefficients. Outside of the cylinder 
(r > r cy i) the magnetic field can be represented as the 
sum of the incident field and the scattered field: 



-Hout — -ffiNC + H s 



(3) 



The incident field under the QS approximation is linear 
with respect to r: 



Hmc = H + h Q r cos(0 - <f> ) 



(4) 



where ho is a coefficient. The constant H does not af- 
fect the results and is for simplicity set to zero in the 
remainder. The scattered field is 



H s (r,<j>) = 



-H 



exp(in</>) 



(5) 



where c„ are coefficients to be determined. 

The boundary conditions across the coating are as fol- 
lows. If the coating is nonmagnetic 2 and thin, the mag- 
netic flux passing through it can be neglected and conse- 
quently the tangential component E<j, of the electric field 
is assumed to be continuous across the coating: 



Substituting 



into (6), we have 



Ea, 



1 dH 
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^in dH out 



dr e out dr 

which links the coefficients a n , c„ as follows: 
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(6) 
(7) 

(8) 

(9) 
(10) 



The jump of the magnetic field is not neglected, as it 
corresponds to a current layer that may be appreciable if 
|e s heii| is large (e.g. for plasmonic materials): 



H out - H in = iui e she ii(<^) Ea, Ar(<j>) 
With E<j, defined by (7), equation (11) becomes 



e s heii(<£) A? ,^ dH iD/out 



a/out 



dr 



(ii) 



(12) 



where the "in/out" subscript in the right hand side en- 
compasses two equally valid expressions: one with e; n and 
Hi n and another one with e ou t and H out - 

From superfocusing at a point to fields on the 
cylinder. Let the desired field distribution around a 
certain focusing point r = rf > r cy \, <f> = <f)f be given 
by (1). This behavior of the field defines the expansion 
coefficients c„: 



H f 

c « = 9n, ±1; c±i = 



1 



h r f exp(=Fi</»o) 

rf l 

(13) 

<pf) exp(in<p)d<p arc the 



where g n = (2ir) 
Fourier coefficients. 

With the coefficients c„ so defined, one can now evalu- 
ate a n from (9), (10), then the jump of the magnetic field 
and consequently, from (12), the required parameters of 
the coating. 



e she ii(<A) Ar(4>) = 



(H out (r cy i) - H m (r cy i)) 
dH in (r cyl )/dr 



(14) 



1 The magnetic field can itself be viewed as a (scaled) stream func- 
tion. 



2 True for all natural materials at optical frequencies. 
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Fig. 2. esheii Ar/r cy i vs. angle for Hf = 2, r f = l.2r cy \ 
and A<f> = 2tt/16, n max = 12 harmonics. 



This result is valid for thin nonmagnetic shells and also 
holds for the full electrodynamic problem, although the 
coefficients a n , c n implicit in (14) are different in the 
wave case. 3 

Note that Eq. (14) defines the product of the per- 
mittivity and thickness of the shell rather than each of 
these quantities separately. This reflects the impedance- 
like nature of the thin-shell boundary condition: Z = 
E^/(H out — H ln ) = i/ (u)e s b e nAr) has the physical mean- 
ing of impedance. 

To summarize, the algorithm proceeds as follows: 

1. Choose r cy i, rj and the maximum number n max of 
cylindrical harmonics. 

2. Compute c n from (13), for \n\ < n max . 

3. Find a n from (9), (10). 

4. Compute H out (r cyi ), H in (r cy i) and dH in (r cyl )/dr 
using the cylindrical harmonic expansion with the 
coefficients a n , c n found previously. 

5. Find ChenW) Ar(0) from (14). 

Changing the desired field variation g(cj>) and the ampli- 
tude Hf, one obtains a very rich variety of solutions for 
the shell parameters (14). To illustrate some of the possi- 
bilities, consider the angular distribution of e s heii Ar/r cy i 
at r — rj — 1.2r cy i- For relatively small values of Hf, 
the sign of e s heii varies (Fig. 2, Hf — 2, A(p = 27r/16, 

^max — 12). 

Higher values of Hf require a strong plasmonic reso- 
nance with e s hoii(</ ) ) consistently negative (Fig. 3, Hf 
10. Umax = 20 and A<fi = 2n/18). Here the deriva- 
tive g'(4>) was chosen as a Gaussian peak g'(<f>) — 
exp(— (f> 2 /A(p 2 ), where parameter A</> controls the width 



3 As in [8], in the wave case small negative values of e"jj ell (</>) in 
(14) are mathematically possible but could in practice ignored. 
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Fig. 3. e sh eiiAr/r cy i vs. angle for Hf = 10, r f = 1.2r cy \ 
and A<p = 27r/18, n max = 20 harmonics. 



of the peak. This angular distribution of the magnetic 
field leads to a respective peak in E r (<fi) (Fig. 4 and 
Fig. 5). The fields were computed semi-analytically with 
(3)-(5) and also numerically using finite-difference (FD) 
analysis on regular polar grids. The semi-analytical so- 
lution in its present form is valid only in the absence of 
losses (e" heU = 0) and, as Fig. 5 shows, is in that case 
close to the numerical solution. (The discrepancy is due 
primarily to the finite thickness of the shell Ar = 0.02r cy i 
in FD simulations, and to the limited number of cylindri- 
cal harmonics in the semi- analytical treatment.) Losses 
affect the amplitude of the peak but not its sharpness 
(Fig. 5, e" holl = 0.1). Convergence of FD results was ver- 
ified by running the simulation on different grids (e.g. 
almost coinciding dashed and dotted lines corresponding 
to two different grids, Fig. 5). 

In summary, "superfocusing" of light by nanoshcll 
scatterers, with the sharpness of the focus in principle 
unlimited, has been demonstrated analytically and nu- 
merically. The shell is designed in two stages: (i) ana- 
lytical continuation of the desired behavior of the field 
at the focus to the boundary of the scatterer; (ii) finding 
the distribution of the dielectric permittivity and thick- 
ness of the shell that would produce the required field 
on the surface. Recently proposed focusing by "super- 
oscillations" [4] also falls into this framework, but with 
the waveform at the focus cleverly chosen to contain only 
traveling wave components. 

The procedure can be generalized to arbitrary shapes 
of the scattering shell and to 3D by still using cylindrical 
/ spherical harmonic expansion in the outside region (as 
in T-matrix methods [15]) and numerical methods (e.g. 
FEM) inside. 

Numerical examples of very sharp focusing presented 
in this Letter should be viewed primarily as proof of con- 
cept; further improvement of the focusing effects should 
definitely be possible with sophisticated numerical op- 



4 



" = 20, r, = l.2r cyl SliperfOCUS 




Fig. 4. Surface plot of E r reveals "super-focuses" inside 
and outside the shell. No losses. Parameters: Hf = 10, 
"max = 20 harmonics, r/ = 1.2r cy i, Atfi — 27r/18, FD 
grid n r xn<f, = 300 x 720. 
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Fig. 5. E r vs. angle at the radius of the focus rf = 
1.2r cy i. Same parameters as in the previous figure. 
Losses (e" ho n = 0.1) reduce the magnitude of the focusing 
peak but do not significantly affect its sharpness. Solid 
line: semi-analytical solution (20 harmonics). Dashed 
and dotted lines on top of one another: no losses, FD 
solutions for grids n r x = 150 x 480 and 300 x 720. 
Empty squares: FD solution, e" = 0.1, grid 150 x 480. 

timization techniques (e.g. adaptive goal-oriented finite 
element analysis [14]) that will treat not only the physi- 
cal properties of the shell but also its geometric shape as 
adjustable parameters. Finally, it would be interesting 
to explore similar ideas for superfocusing of surface plas- 



mon polaritons propagating on surfaces with judiciously 
chosen patterns [4] . 

I thank N. I. Zhcludcv, S. I. Bozhevolnyi and the 
anonymous reviewer for very interesting and helpful com- 
ments. 
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